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Abstract. The acceleration of energetic particle trans- 
port in high amplitude magnetosonic and Alfvenic turbu- 
lence is considered using the Monte Carlo particle simu- 
lations which involve integration of particle equations of 
motion. We derive the momentum diffusion coefficient Dp 
in the presence of anisotropic turbulent wave fields in the 
low-/? plasma, for a flat and a Kolmogorov-type turbulence 
spectrum. We confirm the quasilinear result (cf. Schlick- 
eiser & Miller (1997)) of enhanced values of due to 
transit-time damping resonance interaction in the pres- 
ence of isotropic fast-mode waves as compared to the case 
of slab Alfven waves of the same amplitude. The relation 
of Dp to the turbulence amplitude and anisotropy is in- 
vestigated. 

Key words: cosmic rays - magnetohydrodynamic turbu- 
lence - interstellar medium - Fermi acceleration 



• I— I , 1. Introduction 
V . 

5—1 ' l^s-iiy astronomical objects (extragalactic radio sources, 
\ supernova remnants, solar flares) emit radiation with non- 
thermal spectra. These emissions are often connected with 
the existence of a hot turbulent magnetized plasma pro- 
viding conditions for particle acceleration by MHD turbu- 
lence. R was first shown by Hall and Sturrock (1967) and 
by Kulsrud and Ferrari (1971) that charged particles can 
be accelerated by MHD turbulence having wavelengths 
long compared to the particle gyration radius. For ex- 
ample, the importance of a Fermi-like acceleration mech- 
anism was proved for extragalactic radio sources (Burn 
1975; De Young 1976; Blandford and Rees 1978; Achter- 
berg 1979 and Eilek 1979) and for the second phase ac- 
celeration in solar flares (Melrose 1974; Ramaty 1979). 
The interaction between the waves and particles is de- 
termined primarily by resonances w — = nfi between 



particles characterized by their parallel velocity W|| and gy- 
rofrequency 17 and undamped waves characterized by their 
frequency uj and wavenumber fc El. At such interactions en- 
ergies vary in a diffusive way, and over a long timescale 
a net energy gain results in the process of stochastic ac- 
celeration. For n 7^ 0, a particle in its parallel rest frame 
(w|| = 0) resonates with those waves that it sees at an inte- 
gral multiple of a gyrofrequency and for n = at zero fre- 
quency. Then, a particle feels a net force which is not aver- 
aged out by the phase mixing (Stix 1962). R is a magnetic 
analogy of the Landau damping and is associated with 
the first-order change in (5|B|. This type of coupling is not 
observed in the case of Alfven waves but becomes impor- 
tant for the magnetosonic waves containing compressive 
components of (5B. The interaction between the particle 
magnetic momentum and the parallel gradient of the mag- 
netic field is called transit-time damping (cf. discussions 
by Lee and Volk 1975, Achterberg 1981, Miller et al. 1996). 
Recently, Schlickeiser & Miller (1997) presented a quasi- 
linear derivation of cosmic ray transport coefhcients in the 
presence of MHD waves, including the isotropic fast-mode 
turbulence. For the isotropic Kolmogorov turbulence they 
demonstrated that the Fokker-Planck coefficients depend 
both on the transit-time damping and the gyro-resonance 
interactions. For cosmic ray particles with u 3> Va and a 
vanishing turbulence cross helicity the momentum diffu- 
sion coefficient in the fast-mode turbulence is mainly de- 
termined by the transit-time damping contribution, lead- 
ing to a more effective stochastic acceleration in compar- 
ison to the process in the presence of the pure Alfvenic 
turbulence. 

The aim of the present paper is to study the momen- 
tum diffusion coefficient Dp in the presence of non-linear 
((5_B > Bq) magnetosonic and Alfven waves where the 
quasilinear approximations do not apply. We also study 
the influence of the degree of wave propagation anisotropy 
and of the wave spectral slope on the value of Dp. To 
consider these problems the Monte Carlo simulations in- 
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volving derivations of particles trajectories in the space 

filled with finite amplitude fast-mode and Alfvcn waves 
are applied. Anisotropic wave distributions are modeled 
by choosing their wave vectors from finite opening cones 
directed along the mean magnetic field Bg - For our simula- 
tions we adopt fast-mode or Alfven mode turbulence with 
flat ((7=1) and Kolmogorov-type {q = 5/3) power law 
spectrum cx k~'^ above the minimum wavenumber k^in- 
We confirm a substantial increase of Dp for the fast-mode 
waves in comparison to the Alfven waves of the same am- 
plitude if the nearly perpendicularly propagating waves 
[k _L Bo) are included. 

2. Quasi-linear momentum diffusion coefHcient 

The quasi-linear theory treats the effect of the weakly per- 
turbed magnetic field as perturbations of orbits of parti- 
cles moving in the average background field. Schlickeiser 
(1989) has considered the quasilinear transport and ac- 
celeration parameters for cosmic ray particles interact- 
ing resonantly with Alfven waves propagating parallel to 
the average magnetic field. The transport equation can 
be derived from the Fokker-Planck equation by a well- 
known approximate scheme (Jokipii 1966, Hasselmann & 
Wibberenz 1968) which is commonly referred to as the 
diffusion-convection equation for the pitch-angle averaged 
phase space density. Into the equation the electromag- 
netic fields generated by MHD waves enter through the; 
Lorentz force term. For fast {v ^ Va) cosmic ray parti- 
cles, a vanishing cross helicity state of the Alfven waves 
and the power-law turbulence spectrum 



the momentum diffusion coefficient 
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Recently, Schlickeiser & Miller (1997) considered cosmic 
ray particles interacting with oblique fast-mode waves 
propagating in a low-/? plasm. In the cold plasma limit the 
fast and slow magnetosonic waves merge to the fast-mode 
waves with the same dispersion relation = V^k"^. For 
such waves they demonstrated that the rate of adiabatic 
deceleration vanishes, and the momentum diffusion coeffi- 
cient and spatial diffusion coefficient fc|| can be calculated 
as pitch angle averages of the functions determined by 
the non-vanishing Fokker-Planck coefficient Dju^u and by 
Dpp = e^p^£)pp(/i) , where e = Va/v. Adopting isotropic 
fast-mode turbulence with a power-law turbulence spec- 
trum they obtained 
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which is the sum of transit-time damping (/y) and gy- 
roresonance interaction (/g) contributions. The form of fx 
admits the pitch angle scattering by transit-time damping 
of super- Alfvenic particles with pitch-angles contained in 
the interval e <| /i |< 1. In the interval | |< e, where no 
transit-time damping occurs, the gyroresonance {n ^ 0) 
interactions provide a small but finite contribution to the 
particle scattering rate. As a result the momentum dif- 
fusion coefficient Dp is predominantly determined by the 
transit-time damping contriburtion and the spatial diffu- 
sion co(?ffic,i(^nt is determined only by the gyroresonance 
interactions. Finally for the fiat (1 < g < 2) and the steep 
{q > 2) fast-mode turbulence spectra (eq. 2.1) they ob- 
tained 
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respectively. 



(2.3) 



3. Description of simulations 

The approach applied in the present paper is based on 
numerical Monte Carlo particle simulations. The general 
procedure is simple: test particles are injected at random 
positions into a turbulent magnetized plasma and their 
trajectories are followed by integrating the particle's equa- 
tions of motion. Due to the presence of waves, particles 
move diffusively in space and momentum. By averaging 
over a large number of trajectories one derives the diffu- 
sion coefficients for turbulent wave fields. In the simula- 
tions we consider relativistic particles with v ^ Va and 
use dimcnsionless imits (cf Appendix A): SB = SB /Bo 
for magnetic field perturbations, l/fio for time, k/kres for 
wave vector and Pgflo for the momentum diffusion coeffi- 
cient. 

3.1. The Wave Field Models 

In the modelling we consider a superposition of 384 MHD 
waves propagating oblique to the average magnetic field 
Bo = Bo^z- The wave propagation angle with respect to 
Bo is randomly chosen from a uniform distribution within 
a cone ('wave-cone') along the mean field. For a given 
simulation two symmetric cones are considered centered 
along Bo, with the opening angle 2a, directed parallel 
and anti-parallel to the field direction. The same number 
of waves is selected from each cone in order to model the 
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symmetric wave field. Related to the wave 'i' the magnetic 
field fluctuation vector SB^''^ is given in the form: 

= sin(A;W • r - w^i) . (3.1) 

The electric field fluctuation related to a particular wave 
is given as SE^^^ = — A 6B^'^^ where V^'^^ is a wave 
velocity. 

For Alfven waves (A) we consider the respective dis- 
persion relation 

a;i = fcfyi , (3.2) 

where Va = Bo/\/^JTp is the Alfven velocity in the field 
Bo- The wave magnetic field polarization is defined by the 
formula 

6BA = SBA{k,L0A){kx e;,)k]^^ . (3.3) 

In low-/? plasma the fast-mode magnetosonic waves (M) 
propagate with the Alfven velocity and the respective re- 
lations are: 

= k^Vi (3.4) 

SBM = SBM{KujM){kx{kxe^})k-^k2^ . (3.5) 

For all our simulations we adopt Va = 10~^c. One should 
be aware of the fact that the considered turbulence model 
is unrealistic at large 6B and the present results can not 
be considered as the exact ones. In particular, in the pres- 
ence of a finite amplitude turbulence the magnetic field 
pressure is larger than the mean field pressure and the 
wave phase velocities can be greater than the assumed 
here VaIBo). 

3.2. Spectrum of the turbulence 

In our simulations we consider a power turbulence spec- 
trum where the irregular magnetic field, obtained from the 
energy density which is defined in equation (2.1), in the 
wave vector range {kmin, kmax) can be written 

/ h \ -g/2 

5B{k) = 5B{kmin)(-r ) (3.6) 

where fcrnm = 0.08(A;TOaa; = 8.0) corresponds to the con- 
sidered longest (shortest) wavelength and q is the wave 
spectral index. In the present simulations we consider the 
flat spectrum with q = \ and the Kolmogorov spectrum 
with q = 5/3. We included the flat spectrum because of 
our earlier simulations (Michalek & Ostrowski 1996) when 
we considered the Alfv'en waves with g = 1. On the other 
hand such kind of turbulence spectrum is very convenient 
for numerical simulations due to presence of a big num- 
ber of short waves. For the flat spectrum wave vectors 
are drawn in a random way from the respective ranges: 
2.0 < fc < 8.0 for 'short' waves, 0.4 < fc < 2.0 for 'medium' 
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waves and 0.08 < k < 0.4 for 'long' waves. The respective 
wave amplitude are drawn in a random manner so as to 
keep constant 

384 

[Y^{6Bi^^r]'/'^6B, (3.7) 

i=l 

where SB is a model parameter, and, separately in all 
mentioned wave- vector ranges 

128 256 

[E(JB«)^]^/^.[E(<5S«)^]^/^. 

j=l 1=129 

384 j.„ 

-[E('5^^"n'''-^- (3-8) 

i=257 V3 

Thus the wave energy carried is uniformly divided into the 
all wave- vector ranges. As a second "realistic one" turbu- 
lence model we consider one with the Kolmogorov spec- 
trum. The observed spectra in the interplanetary space 
often take a such form (Jokipii 1971). In that case all 
wave vectors are drawn in a random manner from the 
whole considered range (0.08 < k < 8.0) but the ampli- 
tudes are fltted according to Kolmogorov law (2.1) 
and to keep the formula (3.7). In a such turbulence spec- 
trum most of the energy is carried by 'long' waves. In the 
discussion below we will consider four different turbulence 
fields: i. Alfven waves with a fiat spectrum - AF, ii. Alfven 
waves with a Kolmogorov spectrum - AK, iii. Fast-mode 
waves with the flat spectrum - MF, iv. Fast-mode waves 
with the Kolmogorov spectrum - MK. 

4. Results 

The results of the simiilations are presented at the figures 
1 and 2. In Fig. 1 on the successive panels, the simulated 
momentum diffusion coefficients Dp for the Alfven and 
the magnetosonic turbulence are presented. The derived 
diffusion coefficients Dp are given for different wave-cone 
opening angles and for different turbulence amplitudes. In 
all cases a systematic increase of Dp with the amplitude 
occurs, but the rate of this increase diminishes at larger 
non-linear SB. We suspect the considered effect reflects 
partial particle trapping in high-amplitude turbulence, de- 
creasing randomness of the particle momentum variation. 
There are no significant differences between results ob- 
tained for both types turbulence spectrum, q = I or 5/3. 
For the Alfvenic turbulence an increase of the wave-cone 
opening angle a provides waves with smaller phase veloci- 
ties and leads to decreasing values of Dp. Additionally, the 
values of the wave vector amplitude fcy decrease leading 
to less efficient coupling mediated by the cyclotron reso- 
nance. The trend is independent of the considered turbu- 
lence amplitude. In the case of fast-mode waves a more 
complicated relation is observed. For the flat turbulence 
spectrum a small increase of the angle a leads to a slight 
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Fig. 1. Variation of the momentum diffusion coefficient Dp versus the perturbation amphtudc SB and the wave propagation 
anisotropy for the flat spectrum and the Kolmogorov spectrum. For comparison results for both the Alven turbulence (thin 
lines) and the fast-mode turbulence (thick lines with indicated simulation points) are superimposed at the lowest panels. 



decrease of Dp due to mentioned weakening of cyclotron 
resonance. However, the appearance of more waves propa- 
gating nearly perpendicular to the ordered magnetic field 
inverts this trend by allowing the transit time-damping 
resonance to become more effective. For the Kolmogorov 
turbulence an increase in opening angles is always followed 
by an increase of Dp, but, again reaches the maximum in 
the presence of perpendicular waves. The described non- 



monotonic variations are clearly seen at low wave ampli- 
tudes in Fig. 1, while for large amplitudes the phenomenon 
is less pronounced. Generally, for large amplitudes of the 
magnetosonic waves the value of Dp does not strongly de- 
pend on the opening angle of the waves. The simulation 
errors can be evaluated from comparison of the results for 
a = 0° of Alfven and fast mode waves, which should coin- 
cide. In Fig. 1 only the clarity the results for a = 0°, 40° 



Cosmic Ray Momentum Diffusion 



5 



and 90° are presented, but the ones calculated by us for 
a = 30° and 60° arc consistent with the above description. 

In order to explain this non-monotonic behavior one 
can relate to the quasi-linear derivations of Schlickeiser 
& Miller (1997). For n = the resonance condition for 
the transit-time damping may be written as = Lv/k^\ 
and, with the dispersion relation (4), Va/u|| = fcy/fc. It 
is clear that for Va v particles can effectively inter- 
act with waves at wide range of only if the waves 
with k± » fc|| are present. This fact explains the effective 
transit-time damping interactions for magnetosonic waves 
with isotropic space distribution of waves. Then, the Dp 
increase results mainly from the presence of waves prop- 
agating quasi-perpendicular to the mean magnetic field. 
In our simulations we adopted highly relativistic particles 
with V = 0.99c » Va, so that this effect is the most pro- 
nounced. 

In order to verify this behaviour in more detail we 
considered the interaction of relativistic particles with 
waves propagating in narrow cones with respect to the or- 
dered magnetic field. We performed simulations for cones 
0° ± 5°, 45° ± 5°, 85" ± 5° and results are shown at Fig. 2. 
Only the presence of waves propagating nearly perpendic- 
ular to the mean magnetic field leads only to significant 
increase of Dp. For these waves the transit-time damping 
resonance is responsible for the wave particle coupling, 
while the waves propagating at small angles or parallel to 
the average magnetic field only contribute to the gyrores- 
onance interaction. 

5. Discussion and summary 

We considered momentum diffusion coefficient Dp in the 
presence of oblique Alfvenic and magnetosonic turbulence 
with amplitudes ranging from small ones up to highly 
non-linear ones. The infiuence of the degree of anisotropy 
and of the waves spectrum slope on Dp was also stud- 
ied. Generally, in all cases a systematic increase of Dp 
with the wave amplitude is observed. Dp does not strongly 
depend on the type of turbulence spectrum. We confirm 
a substantial increase of Dp for the fast-mode waves in 
comparison to the Alven waves of the same amplitude 
if the nearly perpendicularly propagating waves are in- 
cluded. That increase is caused by the transit-time damp- 
ing which occurs in the presence magnetosonic of waves 
containing compressive component of SB. In the presence 
of a small amplitude magnetic turbulence a significant in- 
crease of Dp, on a factor of 10, is achieved. This agrees 
well with the estimate of Schlickeiser and Miller (1997) 
that the coefficient Dp for the isotropic Kolmogorov-type 
spectrum in the presence of magnetosonic waves is about 
a factor Iny^ =o(10) larger than the one for the parallel- 
propagating Alfven waves. In the case of our simulations 
for the Kolmogorov spectra this value is nearly constant 
for a wide range of waves' amplitudes and for isotropic 
wave spectrum our results are consistent with the quasilin- 



(5B = 0.6 




20000 40000 60000 80000 100000 
time 

Fig. 2. Examples of the simulated Dp for anisotropic waves 
propagating close to the chosen pitch angle (0°,45°,85°) for 
both types of turbulence spectrum and 6B = 0.6. 

ear estimates. The approximation worsens for the largest 
non-linear amplitudes (SB/Bq) > 1.5) when the consid- 
ered ratio diminishes on about 50%. 

We also found that for large amplitudes of magne- 
tosonic waves the difference between Dp does not strongly 
depend on the opening angle of the waves. When the 
strength of the compression is high enough most particles 
can be reflected before complete penetrating the region 
of the compression they could only be accelerated due to 
a simple second order Fermi mechanism. As the ampli- 
tude decreases, the number of particles that are reflected 
decrease. For very small amplitudes particles can be re- 
flect at the wave when their parallel velocities in the wave 
frame is about zero. That occurs when v\\ = oj/k\^ what is 
condition for the transit-time damping. Hence, that pro- 
cess could be called small-amplitude Fermi acceleration 
(Achterberg 1981). We should remember that the wave 
phase velocity are expected to be larger in realistic turbu- 
lence than the assumed here Va{Bo). Therefore our results 
for large SB should be consider rather as the lower limits. 

Considerations of the MHD waves propagation oblique 
to the mean magnetic field (e.g. Tademaru 1969, Lee and 
Volk 1975) show that such waves are subject to effective 
processes dissipating their energy. Because of that the con- 
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sidered here effects of the momentum diffusion enhance- 
ment due to fast-mode waves can occiir only in a vohimc 
with acting the turbulence generation force. For example, 
in vicinity of the strong shock or in a region of the mag- 
netic field reconnection the required fast-mode perpendic- 
ular waves arc expected to be effectively created. 
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e - particle charge 

^res — 

fc|| - wave vector along mean magnetic field 

p - particle momentum vector 

q - spectral index of waves 

Tg - particle gyro-radius 

Rl =1^1" particle Larmour radius 

V = c^p/e - particle velocity vector 

Va - the Alfven velocity in the field Bq 

W|| - velocity along mean magnetic field 

a openin angle of MHD waves 

7 = (1 — z;^/c^)~^/^ - the Lorentz factor 

i - number of give waves 

k — wave- vector 

K± - transverse (cross-field) diffusion coefficient 

K|l - parallel diffusion coefficient 

TO - particle mass (to = 1 in simulations) 

u> - wave frequency 

fl = eB/'jmc - particle angular velocity 

rj(°) = eBoh-mc 

G - the momentum pitch-angle with respect to Bq 

{fjb = cos O) 

£ - particle energy 



A. Summary of notation 

02 = Dp- momcntiim diffusion coefficients 
B = Bo -|- <5B - magnetic induction vector 
Bo - regular component of the background magnetic 

field 

{Bq = 1 in the simulations) 

5B - turbulent component of the magnetic field 

c - light velocity (c = 1 in simulations) 

Z?^^ Fokker-Planck coefficient 

E electric field vector 

(fc) - magnetic energy density for given wave 
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